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ANALYSES  OF  PLATES  CONSTRUCTED  OF  * 
FIBER-REINFORCED  BIMODULUS  COMPOSITE  MATERIAL 


ABSTRACT 

To  implement  the  structural  application  of  the  recently  introduced  fiber- 
governed  symmetric  compliance  model  for  bimodulus  composite  materials,  both 
classical  closed-form  and  finite-element  solutions  are  being  developed.  This 
paper  summarizes  the  results  obtained  for  deflection  of  single-layer  ortho- 
tropic and  two-layer,  cross-ply  plates  of  the  following  configurations  and 
loadings: 

1.  Thin  elliptic  plates  clamped  on  the  boundary  and  subjected  to 
uniform  pressure. 

2.  Moderately  thick  rectangular  plate  freely  supported  on  the  boundary 
and  subjected  to  sinusoidally  distributed  pressure. 
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NOMENCLATURE 

stretching  stiffness  (i,j=l,2,6) 
plate  semi-axes 

bendir.g-stretching  coupling  stiffness  ( i , j=l  ,2,6) 
bending  stiffness  ( i , j=l  ,2,6) 
a(  )/ax 

finite-element  force  components  (i=l  ,2 , . . . ,5) 
total  plate  thickness 

stiffness  coefficients  in  the  finite-element  formulation 
linear  differential  operators  defined  in  equations  (2) 
stress  couples  and  stress  resultants  (i=l,2,6) 
finite-element  interpolation  functions 
number  of  nodes  per  element,  number  of  layers  in  laminate 
normal  pressure,  intensity  of  normal  pressure 
thickness  shear  stress  resultants 


This  paper  is  to  be  presented  at  the  Symposium  on  Mechanics  of  Bimodulus 
Materials,  sponsored  by  the  Applied  Mechanics  Division,  at  the  ASME  Winter 
Annual  Meeting,  New  York,  NY,  Dec.  1979  and  will  be  published  in  the 
symposium  volume. 
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plane-stress  reduced  elastic  stiffness  (1  ,j  = 1 ,2,6) 
for  tension  ('<*1)  or  compression  (k=2) 
for  layers  1 or  2 ( 1=1  or  2) 
elastic  compliance  coefficient  (i,j=l,2,6) 

finite -element  matrix  coefficients 
midplane  displacements  in  x,y,z  directions 
total  displacements  in  x,y  directions 
coefficients  in  expressions  for  u,v,w 
in-plane  rectangular  position  coordinates 

znx/h’  zny/h 

thickness-direction  position  coordinate,  measured  from 
midplane 

2 coordinates  of  neutral  surfaces  based  on  e and  e 

x y 

fiber-direction  strain 
strain  (j=l  ,2,6) 
stress  (i=l  ,2,6) 

angle  between  fiber  direction  of  an  individual  layer  and 
reference  direction  for  the  laminate 

slope  functions 

middle  plane  of  the  plate 

a typical  finite  element 


INTRODUCTION 


An  elastic  plate  is  an  important  structural  component  in  a wide  variety 
of  engineering  systems.  Thus,  it  is  not  surprising  that  the  first  modern 
development  of  the  basic  constitutive  equations  of  bimodulus1  materials  by 
Ambartsumyan  (IJ2  was  followed  only  one  year  later  by  Shapiro's  analysis  (2j  of 
static  deflection  of  a circular  plate  constructed  of  such  a material  and  sub- 
jected to  a pure  radial  bending  moment,  however,  in  his  analysis,  he  used 
Love's  stress-function  formulation  rather  than  plate  theory. 

Unfortunately,  there  have  been  relatively  few  more  recent  analyses  of 
bending  of  bimodulus  plates.  Notable  exceptions  are  the  series  of  papers  by 
Kamiya  (2*1)-  ln  (2).  he  treated  large  deflections  (geometric  nonlinearity) 
of  circular  elates,  using  an  iterative  finite-difference  technique,  while  in 
(4)  , he  applied  the  energy  method  to  large  deflections  of  a rectangular  plate. 
The  effect  of  thickness  shear  deformation  on  the  linear  problem  was  treated 
in  (5). 

Ponomarev  (6)  considered  bending  of  a square  plate  made  of  a nonlinear 
elastic  material  of  a slightly  more  general  nature  than  bilinear,  namely  one 
with  a third-degree  polynomial  stress-strain  relation  such  that  the  ratio  of 
tension  stress  to  compression  stress  at  the  same  absolute  value  of  strain 
remains  constant. 


All  of  the  analyses  mentioned  above  are  limited  to  isotropic  materials 
with  different  properties  in  tension  and  compression.  Apparently,  the  first 


1 Here  a bimodulus  material  is  understood  to  be  one  having  different  elastic 
propertiesln  tension  (T)  and  compression  (C).  Materials  with  same  properties 
in  T and  C are  referred  to  herein  as  ordinary  materials. 

2 Underlined  numbers  in  parentheses  designate  References  at  end  of  paper. 
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analysis  considering  bimodulus  materials  other  than  isotropic  is  due  to  Jones 
and  Morgan  (7),  who  presented  a closed-form  solution  for  cylindrical  bending  of 
a thin,  cross’-nly  laminate  of  finite  width  but  infinite  length.  They  presented 
numerical  results  for  examples  of  both  antisymmetric  and  more  generally  un- 
Symmetric  cross-ply  constructions. 

Apparently,  the  first  analyses  of  anisotropic  bimodulus  plates  finite  in 
both  directions  are  those  recently  carried  out  by  the  present  authors  and  their 
colleagues  and  summarized  here.  The  planform  geometries,  boundary  conditions, 
loadings,  lamination  schemes,  degree  of  nonlinearity,  and  methods  of  solution 
are  listed  in  Table  1. 


Table  1.  Specific  Plate  3ending  Problems  Considered 
For  Anisotropic  Bi.modulus  Materials* 


[Problem 

iNo. 

Planform 

Geometry 

Boundary 

Conditions 

Lamination 

Scheme 

Degree  of 
Nonlinearity 

Method  of 
Solution 

Elliptic 

thin 

Clamped 

Unidirectional 

orthotropic 

Linear 

1 

Closed  form  (8) 

* 

Elliptic 

thin 

Clamped 

Unidirectional 
orthotropic  & 
cross-ply 

von  Karman 
geometric 
nonl ineari ty 

Perturbation 

(9) 

1 3 

Elliptic 

thin 

Clamped 

Arbitrary 

cross-ply 

Linear 

Closed  form 

4 

Rectangular 

thick 

Freely 

supported 

Unidirectional 
orthotropic  & 
cross-ply 

Linear 

Closed  form 
and  finite 
element 

5 

Rectangular 

thick 

Freely 

supported 

Cross-ply 

Linear 

Finite  element 

GOVERNING  DIFFERENTIAL  EQUATIONS 

Classical  small-deflection  theory  of  thin,  unsymmetrical ly  laminated,  ani- 
sotropic plates  was  originated  by  Reissner  and  Stavsky  (K)).  Apparently,  the 
first  display  of  such  a theory  for  the  completely  arbitrary  anisotropic  lamin- 
ated plate  was  due  to  Whitney  and  Leissa  These  same  equations  govern  the 

small  deflections  of  thin,  laminated,  bimodulus  anisotropic  plates  and  thus 
they  are  presented  here  in  concise  form  for  completeness: 


Ljl  1-12  1-13 

V 

’ o ' 

L 1 2 1-22  1-23 

V 

. = ■ 

0 

_L13  1-23  1-33  . 

>w- 

p , 

where  the  L „ are  symmetric  linear  differential  operators  (L„  = L „1  defined 

as  follows:02  2a  aB 


l12  = A16cij;  + (Aj2  + AG6)dxdy  + A26d* 


♦ 

The  program  can  handle  arbitrary  lamination  arrangement  and  boundary  cond- 
itions. However,  the  results  for  only  cross-ply  laminates  and  freely  supported 
edges  are  included  here. 

* Uniform  loading  is  considered  in  all  .ases,  except  in  Problem  4,  in  which 
sinusoidal  loading  is  used. 
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Ljj  2 - qm^  * 3B16dj|dy  - (B12+  230fl)Jxd^  - B2(idy 

L2;  = A60d2  + 2A2Cdxdy  + A22d2  (2) 

L23  2 - * (8,2  + 2B6c)d^dy  - 332bdxd2  - B22d3 

L33  = 0ud^4016d3dy  + 2(D12  + 2066)d=dy  + «D26dxd3  + D22d£ 


Here  u,  v,  w are  the  midplane  displacements  in  the  x,  y (in-plane)  and  z 
(normal)  directions,  dx  = d(  )/dx,  p is  the  normal  pressure,  and  the  plate 
stiffnesses  are  defined  as 


rh/2 

(1 ,z,:2)q1j(z)dz 

-h/2 


(1.3-1 ,2.6) 


(3) 


Here  h is  the  total  thickness  of  the  laminate  and  z is  measured  from  the  mid- 
plane by  convention.  The  Q..(z)  are  the  plane-stress  reduced  elastic  stiff- 
nesses defined  as  follows:  J 


(o.)  = CQ1J3CeJ)  (1,3-1 .2.6)  (4) 

Here  the  contracted  notation  of  composite-material  theory  (1_2)  is  used.  Thus, 
subscripts  1 and  2 refer  to  normal  stress  (or  strain)  action  in  the  x and  y 
directions  and  6 refers  to  shearing  stress  (or  strain)  action  with  respect  to 
the  x,y  axes.  Due  to  thermodynamic  considerations,  the  [Q,.]  matrix  is 
symmetric. 

PLATE  STIFFNESSES  FOR  BIMODULUS-MATERIAL  LAMINATES 

Even  in  laminates  made  of  ordinary  materials  (those  having  the  same 
elastic  properties  in  tension  and  compression),  the  stiffnesses  Q-ji  are  piece- 
wise-constant  functions  of  the  thickness-direction  coordinate  z.  the  indivi- 
dual layers  generally  consist  of  unidirectional  fiber-reinforced  composite 
material  which  is  orthotropic  with  respect  to  its  material -symmetry  axes:  the 
fiber  direction  and  the  two  directions  orthogonal  to  it.  Thus,  when  the 
individual  layers  are  oriented  parallel  to  either  the  x or  y axes  of  the  plate, 
all  of  the  A, j , Bjj,  and  Dy  stiffnesses  involving  shear-normal  coupling  (all 
of  those  with  subscripts  16  and  26)  identically  vanish. 

Furthermore,  when  a laminate  is  composed  of  individual  layers  of  ordinary 
material  arranged  symmetrically  with  respect  to  its  midplane,  al 1 of  the  so- 
called  bending-stretching  coupling  stiffnesses  ( 3 i -i ) vanish.  Thus,  a single- 
layer ordinary  composite  has  no  By  terms  present/  However,  as  discussed  in 
(7)  and  (8),  even  a single-layer  bimodulus  composite,  as  well  as  a so-called 
symmetrically  laminated  multi-layer  laminate,  has  certain  B-j j terms  present  by 
virtue  of  the  nature  of  the  bimodulus  action.  Thus,  the  terms  symmetric,  anti- 
symmetric, and  unsymmetric  have  no  significance  for  bimodulus  laminates. 

A popular  lamination  scheme  is  the  so-called  balanced  angle-ply  laminate. 
This  is  one  consisting  of  an  even  number  of  layers  of  identical  material  and 
thickness  and  having  an  equal  number  of  layers  oriented  at  an  angle  +6  and  -9 
with  respect  to  an  arbitrary  reference  direction..  When  such  a laminate  is 
made  of  ordinary  material,  the  shear-normal  coupling  effects  (Q16  and  Q26 
terms)  for  +9  are  exactly  balanced  by  the  ones  at  -9.  Hence,  A16  = A26  = 0 
and  the  term  "balanced".  In  contrast,  if  this  same  laminate  is  made  of  bi- 
modulus material  and  undergoes  sufficient  bending  action  that  the  fiber-dir- 
ection strains  at  the  top  and  bottom  surfaces  of  the  laminate  are  opposite  in 
sign,  then  the  terms  at  +9  are  not  balanced  by  those  at  -9.  Thus,  for  this 
bimodulus  laminate,  AI6  and  A26  do  not  vanish  and  the  term  balanced  is 
inaporopriate. 

Even  in  the  case  of  a so-called  symmetric  balanced  angle-ply  (SBAP) 


-6- 

laminate,  such  asone  havinci  a lamination  scheme  0/-0/-9/0,  of  ordinary  composite 
material,  even  though  Aib  and  A,,,  vanish,  D1;,  and  D_,0  do  not.  A detailed  de- 
rivation of  general  expressions  for  t),(>  and  02b  for  such  a laminate  with  an 
arbitrary  even  number  of  layers  is  given  in  (i_3).  !n  view  of  tiiese  considera- 
tions, it  is  clear  that  for  an  angle-ply  laminate  of  bimodulus  material,  there 
are  no  vanishing  plate  stiffnesses. 

A listing  of  all  of  the  plate  stiffnesses  that  vanish  for  a variety  of 
lamination  arrangements  of  both  ordinary  and  bimodulus  composite  materials  is 
presented  in  Table  2. 

Table  2.  Plate  Stiffnesses  for  Some  Popular  Lamination 

Arrangements  of  Ordinary  and  3imodu!us  Composite  Materials 


1 

Name  of 

Example 

' List  of  Vanishing 

Plate  Stiffnesses  ' 

Case 

Arranoement 

Laminate 

| Ordinary  Material 

' Si  modulus  Material. 

1 

j Aligned  single 

1 ply  or  al igned 

1 parallel  ply 

0 or 

90/90 

A)6 ,A26.Dig»32S* 

| all  Blj 

Al6*A26.B16*  1 

B26  » B 1 6 »B2 6 

| 

; 

| Cross-ply: 

| 

! 2 

1 

1 n = odd 

0/90/0 

I 

Al6*A26>Dl6*^26« 
all  Bij 

1 a16,a26,b16,  | 

j B26>°16*D26 

! 3 1 

— 

i n = even 

0/90 

a16,a26,b16, 

B26iB16»B26 

A16*A26>B16*  1 
B2S»D16*D26 

j 

| 

Alternating 

angle-ply: 

1 

1 

I | 

| 4 

n = odd 

30/-30/30 

A16,A26.  al1  S.  ■ 

None 

| 5 

n = even 

10/-10 

A16»A26»B11»B12» 
B22  »^1 6 »B26 

None 

1 

1 6 

SBAP 

30/-30/-30/30 

A16,A26.  all  B.j 

None 

i 7 1 

Quasi-lsotropic  ; 

I 

0/-45/+45/90 

A16 ,A26 ,DI6  ,D26 

None 

j 8 

Q/+  8 

-30/0/30 

[ A16,A26,Bl j ,B12, 
B22  >°16  »°26 

None  ■ 

9 

(0/+  e)s 

-30/30/0/30/-30 

A;6,A26,  all  B^j 

None 

The  remainder  of  this  discussion  is  limited  to  the  particular  bimodulus- 
material  model  originated  by  Bert  (]_4_, T_5_)  and  known  as  tne  symmetric  fiber- 
governed-ccmpl i ance  model.  This  model  can  be  expressed  mathematically  as 
follows:  , , 

Uj)  = CSJ1]to1>  (5) 

where  the  compliance  takes  on  different  sets  of  values  as  follows: 


[Sji] 


fv“’ 

if  Of>0 

«V(C> 

if  Oj  < 0 

(6) 


Here  af  denotes  the  fiber  stress  along  its  axis. 


The  symbol  n denotes  the  total  number  of  layers  in  the  laminate. 
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ti 


To  implement  the  material  model  in  plate  theory,  we  introduce  a third  sub- 
script ( k ) to  the  0, , plane-stress  stiffness  matrix  in  which  k-1,2  denote  ten- 
sion and  compression'pronerties,  respectively.  Tliett  by  invoking  the  well-known 


respective  individual  layer, 

one 

obtains  the 

appropriate  Q^: 

( 

'ro  „]  if 

e - 

> 0 

[Q  ] = 

r iji" 

f 

L ijkj 

1 

KjA  if 

Ef 

< 0 

(7) 


where  c{  denotes  the  fiber-di rection  strain  at  any  arbitrary  point.  The  locus 
of  points  at  which  ef  * 0 is  tradi tional ly  called  the  neutral  surface. 


Dependin']  upontheplate  geometry,  boundary  conditions,  and  loading,  even 
in  the  case  of  a single-layer  plate,  the  neutral  surface  in  general  is  not  a 
horizontal  p’ane  surface.  Since  the  properties  (Qj,.)  depend  upon  the  sign  of 
£f  and  the  plate  stiffnesses  (A,-  j ,Ru  ,Dij ) upon  theJQf1,.  as  in  equations  (3), 
it  is  clear  that  in  general  a bi modulus  plate  is  nonhor.ogeneous  in  its  plane, 
i.e.,  the  plate  stiffnesses  depend  upon  position  (x,y).  It  is  noted  that, 
unlike  a plate  with  linearly  tapering  thickness,  here  the  plate  stiffness  in 
general  is  not  only  not  a simple  function  of  position,  its  functional  form  is 
not  even  known  a priori.  In  this  regard,  there  is  a qualitative  analogy  be- 
tween a bimodulus  problem  and  an  elastoplastic  problem,  for  which  the  elasto- 
plastic  boundary  is  unknown  a priori. 


From  the  above  discussion,  it  is  clear  that  in  order  for  a finite  element 
to  be  applicable  to  any  arbitrary  combination  of  planform  geometry,  boundary 
conditions,  and  loading  for  a bimodulus  laminated  plate,  it  must  have  a full 
array  of  plate  stiffnesses  (none  zero). 

SMALL  DEFLECTIONS  OF  THIN,  CLAMPED  ELLIPTIC  PLATES 


In  studying  structural  mechanics  problems  involving  new  classes  of  mat- 
erials, it  is  often  quite  instructive  to  obtain  closed-form  solutions  for  cer- 
tain special  cases.  Then  there  are  no  questions  regarding  numerical  approxi- 
mations, convergence,  etc.  In  the  case  of  plates  laminated  of  ordinary  mat- 
erials, the  closed-form  solutions  due  to  Kicher  (J_6)  for  elliptic  plates  and  to 
Whitney  and  Leissa  (1_1_)  for  rectangular  plates  are  most  outstanding.  Thus,  it 
was  decided  to  investigate  the  applicability  of  these  forms  of  solutions  to 
plates  laminated  of  bimodulus  materials. 

Kicher’s  solution  (16_)  is  the  closed-form  solution  for  a uniformly-loaded, 
clamped-edge  elliptic  plate  of  cross-ply  construction  with  an  even  number  of 
layers  (Case  3 in  Table  2).  The  form  of  this  solution  is 

u * u [1  - (x/a)2  - (y/b)2](x/a)  ; v = v [1  - (x/a)2  - (y/b)2](y/b) 

u u (8) 

W * w0[l  - (x/a)2  - (y/b)2]2 

Here  the  displacement  coefficients  are  readily  determined  by  direct  substitution 
into  governing  equations  (1). 

In  view  of  the  Kirchhoff  hypothesis,  upon  which  the  present  theory  is 
based,  the  x-direction  normal  strain  at  any  arbitrary  location  (x,y,z)  is 

e,  - u,x  - zw,xx  (9) 

Thus,  If  one  sets  e5  = 0 and  solves  f r z,  one  obtains  the  neutral  surface 
position  znx  associated  with  the  normal  strain  in  the  x direction.  It  can 
readily  be  shown  that  when  this  procedure  is  applied  to  the  Kicher  displacements, 
equations  (8),  the  resulting  expression  for  znx  is  independent  of  x and  y. 
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1 


This  moans  that  for  this  particular  problem,  the  plate  stiffnesses  are  uniform 
throughout  the  plate  surface  of  a single-layer  bimodulus  plate  Provided  that 
the  fibers  are  aligned  in  the  x direction.  (A  similar  conclusion  can  be  reached 
for  z|ly  and  t lius , a single-layer  bimodulus  plate  with  the  fibers  in  the  y 
direction. ) 

For  the  single-layer  plate,  it  can  be  shown  (8)  that  the  plate  stiffnesses 
depend  upon  Z(.zn/h)  in  the  following  manner: 

■ 0/2)0!,,,  * 0,j2>  * W,J2  - 

■ 0/8)01, j,  - 0,J2)  * (Q,j2  • 8,j, H/2/2)  OO) 

V’  * 0/20(0, j,  • 0,J2)  * (0iJ2  - 0,,,)(ZV3) 

The  only  thing  remaining  is  the  determination  of  the  neutral-surface  loca- 
tion Z.  In  principle,  an  equation  in  2 can  be  obtained  by  combining  equations 
(10)  with  the  expression  for  u<j/w0  as  required  by  solution  of  equations  (1). 
Although  the  resulting  expression  is  only  cubic  in  Z,  the  algebraic  structure 
of  the  expressions  for  the  coefficients  of  the  cubic  is  so  lengthy  that  it  was 
found  to  be  more  expedient  computationally  to  obtain  Z by  direct  iteration. 

From  physical  considerations,  only  one  unique  solution  for  Z is  meaningful 
(i.e.  real,  with  ;Z!<  1/2).  In  (8),  numerical  results  were  presented  for  a 
composite  material  with  considerable  bimodulus  action,  namely  aramid-cord/ 
rubber  which  has  the  properties  listed  in  Table  3,  taken  from  experimental  data 
of  Patel  et  al.  ( T7 ) . The  exact  bimodulus  results  were  compared  with  those 
obtained  by  ordinary  theory  using  average  properties.  The  necessity  of  using 
bimodulus  theory  is  clearly  demonstrated  as  shown  in  Figures  1 and  2. 


Table  3.  Bimodulus  Properties  of  Aramid  Tire-Cord/Rubber 
Composite,  Reduced  from  Experimental  Data  in  (17) 


Property 

1 Units 

Tension 

Comoression  1 

Major  Young's  modulus 

| MPa 

3,580 

12.0  | 

i Transverse  Young's  modulus 

MPa 

9.09 

12.0 

Inplane  shear  modulus 

MPa 

3.70 

3.70  1 

(3,83) 

(3.83)  ' 

j Major  in-plane  Poisson's 

ratio: 

0.416 

0.205  1 

| Minor  in-plane  Poisson's 

ratio 

0.01105 

0.205  j 

In  view  of  the  nature  of  Cases  2 and  3 in  Table  2,  one  would  expect  the 
Kicher  displacement  functions  to  be  applicable  to  an  arbitrary  cross-ply  plate 
of  bimodulus  material.  In  ( 1_8) , this  is  shown  to  be  true  and  it  is  implemented 
for  the  "most  unsymmetric"  layer,  namely  the  case  of  n=2.  In  this  case, 
equations  (10)  are  replaced  by 


Vh  = (V2)(QijlI+  Q-i j 22 ) + (Qij21  - QijllJZ,  + <Q,j22  - Qij12)Zy 
B^j/h2  » (l/8)(Q.jn  - Oij22)  + (Qij21  - Qijn)(Z2/2)  ♦ (Qij22  - Q1jI2)(Z2/2) 
D1J/h»  • (1/24) (Q1  j 1 1 + Oij22)  ♦ (Q1J2,  - 01jn)(Zj/3)  + (Qij22  - Qiji2){Z^/3) 


01) 


Here  the  last  subscript  l in  stiffnesses  Qi/ki  denotes  the  layer  number  (either 
layer  1 or  2)  and  the  ijk  subscripts  have  tile  same  meaning  as  before. 


In  the  case  of  cross-oly  laminates,  the  elastic 
layers  having  the  fibers  oriented  in  the  x direction 
direction  neutral  surface,  while  those  for  the  other 
the  y-direction  neutral  surface.  Thus,  to  a certain 


properties  used  for  the 
are  determined  by  the  x- 
layers  are  determined  by 
extent  a two-layer  cross- 
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ply  bimodulus  laminate  is  like  a four-layer,  cross-ply  ordinary  laminate  (with 
undetermined  individual  layer  thicknesses).  Thus,  there  arises  a question  as 
to  the  proper  order  of  the  tensile  and  compressive  reoions.  To  determine  the 
correct  order,  the  following  criteria  were  used  in  (IS): 

1.  For  bendinn-oroducing  loadings,  the  neutral  surfaces  to  which  the 
iteration  procedure  leads  must  lie  within  the  plate  thickness,  i.e. 
iZyj  < 1/2. 

2.  Thd  maximum  plate  deflection  must  be  in  the  same  direction  as  the 
normal -pressure  loading. 

In  the  case  of  an  aramid/rubber  plate,  it  was  found  that  an  order  of 
stress-state  tension  T or  compression  C (listed  from  toD  to  bottom)  of  C/T/T/C 
resulted  in  a negative  maximum  deflection  (uDward  for  a downward  pressure), 
while  an  order  of  C/T/C/T  met  both  of  the  criteria  listed  above  and  thus  was 
judged  to  be  the  correct  solution.  On  the  other  hand,  for  the  case  of  a poly- 
ester/rubber plate,  it  was  found  that  C/T/C/T  resulted  in  a neutral  surface 
lying  outside  of  the  plate,  while  C/T/T/C  gave  results  meeting  both  of  the 
above  criteria. 

In  attempting  to  find  a closed-form  solution  for  an  angle-ply  bimodulus 
plate,  the  second  author  found  a closed-form  solution  for  an  arbitrarily  lamin- 
ated ordinary  laminate  (]_9).  However,  it  did  not  result  in  a constant  value  of 
neutral  surface  and  thus,  apparently  it  cannot  be  extended  to  a closed-form 
solution  of  the  angle-ply  bimodulus  plate.  Alternate  approaches  to  this  prob- 
lem are  currently  being  investigated. 

NONLINEAR  DEFLECTIONS  OF  THIN,  CLAMPED  ELLIPTIC  PLATES 

The  geometrically  nonlinear  midplane  strain-displacement  relations  origin- 
ated by  von  Karman  were  incorporated  into  thin,  laminated,  anisotropic  plate 
theory  by  Whitney  and  Leissa  (11),  although  they  did  not  solve  the  resulting 
nonlinear  partial  differential  equations.  Although  the  Ritz-Galerkin  and 
Rayleigh-Ritz  (energy)  methods  are  probably  equally  popular,  in  (9)  it  was 
elected  to  use  a perturbation  procedure  first  used'for  plates  by  Nash,  Cooley  (20) 
for  clamped,  elliptic  isotropic  plates  under  uniform  pressure.  The  only  change 
in  the  procedure  is  the  additional  iteration  procedure  (necessary  to  take  into 
account  the  bimodulus  stiffnesses)  and  the  additional  elastic  constants  due 
the  bimodulus  action. 

SMALL  DEFLECTIONS  OF  RECTANGULAR  PLATES  INCLUDING 
THICKNESS-SHEAR  DEFORMATION 

It  has  been  demonstrated  (21_)  that  fiber-reinforced  composite  materials 
exhibit  much  larger  thickness  shear  effects  than  do  plates  having  the  same 
geometry  but  constructed  of  homogeneous,  isotropic  materials.  The  explanation 
for  this  is  due  to  the  relatively  low  thickness  shear  moduli  relative  to  the 
in-plane  moduli.  The  numerous  theories  including  these  effects  for  laminated 
anisotropic  plates  were  reviewed  in  (_1_2 ) . It  suffices  here  to  mention  that  the 
two  most  widely  used  of  these  theories  are  those  due  to  Vang  et  al . (22)  and 
Whitney  and  Pagano  (23).  It  was  shown  in  (24)  that  the  latter  theory~Ts  more 
accurate. 

In  [25),  the  modal  shapes  used  by  [26]  for  free  vibration  and  buckling 
analyses  of  shear  deformable  cross-ply  rectangular  plates  are  used  to  analyze 
the  static  deflection  of  similar  plates  constructed  of  bimodulus  composite 
materials  and  subjected  to  a sinusoidally  distributed  normal  pressure.  Again, 
as  in  the  case  of  the  elliptic  plate  problem  discussed  in  the  preceding 
section,  the  two  criteria  for  constancy  of  the  neutral -surface  position  are 
satisfied,  and  thus  the  solutions  are  exact,  closed-form  solutions. 
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FINITE-ELEMENT  FORMULATION  FOR  PLATES  INCLUDING 
TRANSVERSE  SHEAR  DEFORMATION 

While  considerable  effort  has  been  expended  in  the  fini te-element  analysis 
of  isotropic  plates,  only  limited  investigations  of  laminated  anisotropic 
plates  can  be  found  in  the  literature.  Pryor  and  Barker  (2_7) , and  Barker,  Lin 
and  Dara  (2m)  used  the  conventional  displacement  finite-element  method  to  ana- 
lyze thick  laminated  plates.  The  element  has  seven  decrees  of  freedom  (three 
displacements,  two  rotations,  and  two  shear  rotations)  per  node.  Exploiting 
the  symmetries  exhibited  by  anisotropic  plates,  Noor  and  Mathers  (29-31 ) studied 
the  effects  of  shear  deformation  and  anisotropy  on  the  accuracy  and  convergence 
of  several  shear-flexible  displacement  finite-element  models  based  on  a form  of 
Reissner's  plate  theory.  The  analysis  was  limited  to  symmetrically  laminated 
cross-ply  plates  and  the  element  used  involved  80  degrees  of  freedom  per  element. 
The  conventional  finite  element,  when  applied  to  relatively  thick  laminated 
plates,  either  has  failed  to  predict  accurately  the  local  deformations  and 
stresses  of  a plate  under  bending  or  is  too  expensive  to  use  due  to  large  de- 
grees of  freedom  involved  for  even  relatively  simple  problems.  Mau  and  Witmer 
(32) , and  Mau,  Tong,  and  Pian  [22)  have  employed  the  so-called  hybrid-stress 
finite-element  method  to  analyze  composite  elates  including  shear  deformation. 
The  hybrid  elements  have  proven  (see  Gallagher  (3d))  to  have  some  convergence 
problems,  and  in  some  cases  they  give  erroneous  results.  Most  recently,  Panda 
and  Natarajan  ( 3_5 ) used,  following  Mawenya  and  Davies  (36_) , the  quadratic  shell 
element  of  Ahmad,  Irons  and  Zienkiewicz  (37_)  with  the  same  normal  rotation 
through  the  thickness  to  claim  improved  accuracy  over  Mawenya  and  Davies  ( 36 ) . 
The  'thickness  concept*  mentioned  in  there  is  essentially  the  same  as  that  used 
in  the  YNS  theory  (21_) . The  authors  were  primarily  concerned  with  the  accuracy 
of  the  element,  and  no  attempt  was  mace  to  solve  new  problems  for  which  there 
do  not  exist  any  closed-form,  or  exact  solutions.  In  an  effort  to  relax  the 
continuity  requirements  on  the  shape  functions,  Reddy  (38)  recently  devised  a 
simple  finite  element  and  successfully  tested  the  closed-form  solutions  of 
Bert  and  Chen  (39).  Only  investigation  that  concerns  with  the  finite  element 
formulations  of  bimodulus  materials  is  due  to  Crose,  et  al.  (aQ_) . Here, 
following  (38),  we  present  the  finite-element  formulation. 

Consider  a plate  of  constant  thickness  h composed  of  thin  anisotropic 

layers  oriented  at  angles  9i,02,  The  origin  of  the  coordinate  system  is 

located  within  the  middle  plane  (x-y)  with  the  z-axis  being  normal  to  the  mid- 
plane. The  material  of  each  layer  is  assumed  to  possess  a plane  of  elastic 
symmetry  parallel  to  the  xy-plane.  We  shall  denote  the  middle  plane  with  n. 

The  YNS  theory  is  based  on  the  following  assumed  displacement  field: 
u = u(x,y)  + zdx(x,y) 

v = v(x,y)  + zoy(x,y)  (12) 

w = w(x,y) 

where  u,  v,  and  w are  the  displacement  components  in  the  x,y  and  z-directions , 

respectively,  u and  v are  the  in-plane  (stretching)  displacements  of  the  middle 

plane,  and  and  o are  the  slope  functions, 
x y 

The  equations  of  motion  associated  with  the  YNS  theory  are, 
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where  p = p(x,y)  is  the  transversely  distributed  load,  and  N . , M. , Q^,  and  Q 
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A thin  orthotropic  material  with  provision  for  a shift  of  the  neutral  sur- 
face (due  to  different  properties  in  tension  and  compression)  has  the  following 
plate  constitutive  relations: 
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Here  denote  the  stiffness  coefficients  in  the  plate  coordinates  of  the 

2-th  layer  in  tension  (k=l)  or  compression  ( k=2 ) , and  zn  is  the  distance  from 
the  midsurface  to  the  neutral  plane  (which  is  unknown  a priori). 

We  assume,  over  each  element  ne,  the  same  kind  of  interpolation  for  all  of 
the  variables, 

ue  = ? u®N°  , ve  = I v®N*r  , etc.  (n=nodes  per  element)  (18) 
i 1 1 1 * 1 
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where  are  the  element  interpolation  (or  shape)  functions,  ant*  u'(,  and  vV  are 
the  nodal  values  of  u1;  and  vj{  , respectively.  Substituting  (Id)  into  the  first 
variation  of  the  tota.  potential  energy  associated  w'th  (12)  and  collecting  the 


coefficients  of  the  variations  (iu.,  tv.,  etc.)  we  obtain 


[<c]!ac)  = (Fe) 


(19) 


where  Ue}  = {{u  },  (v  },  (w  i,  {«.“},  { d *: ) ) . and  are  given  by 

c c c x y i j 
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(21) 


The  element  stiffness  matrices  are  assembled  in  the  usual  manner,  and 
boundary  conditions  of  the  problem  are  imposed  before  solving  for  (a).  In  the 
present  study,  the  four-node  quadri lateral  element  of  the  serendepity  family 
is  used.  The  element  stiffness  matrices  for  this  element  is  of  order  20x20. 

Figure  3 shows  the  influence  of  the  aspect  ratio  (b/a)  and  side-to-thick- 
ness  ratio  (a/h)  on  the  location  of  neutral  surfaces  for  a single-layer,  iso- 
tropic, bimodulus,  simply  supported  rectangular  plate  subjected  to  sinusoidal 
loading, 

p = PQ  sin^x/a)  sin(ny/b) 

The  following  elastic  properties  are  used: 
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Fig.  3 Neutral -surface  location  vs.  plate  aspect  ratio, 
and  side-to-thickness  ratio  for  single-layered 
rectangular  plate  under  sinusoidal  loading. 
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Fig.  4 Neutral-surface  location  vs.  plate  aspect  ratio 
for  two-layer,  cross-ply  (0°/90°)  square  plate 
under  sinusoidal  loading. 

Note  that  for  b/a  = 1,  the  neutral  surfaces  associated  with  x-  and  y-directions 
coincide  (i.e.,  zpx  = zny). 

Similar  results  are  presented  in  Figures  4 and  5 for  a two-layer,  cross-ply 
(0  790°',  rectangular  plate  under  sinusoidal  loading.  The  bimodulus  properties 
used  are  the  same  as  those  listed  in  Table  3.  Note  from  Fig.  4 that  the 
neutral -surface  location,  znx  is  virtually  unchanged  for  aspect  ratio  greater 
than  1,  while  the  neutral-surface  location,  zny  increases  proportional  to  the 
aspect  raio.  It  Should  also  be  noted  that  the  neutral  surfaces  do  not  coincide 
in  the  cross-ply  case  for  b/a  = 1. 

Figure  6 shows  the  influence  of  the  aspect  ratio,  and  side-to-thickness 
ratio  on  the  transverse  deflection  for  single-layer,  and  two-layer  cross-ply 
problems  d’scussed  above.  The  effect  of  thickness  on  the  deflection  is  more 
pronounced  than  the  efrect  of  the  aspect  ratio. 
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Fig.  5 Neutral -surface  location  vs.  side-to-thickness 
ratio  for  two-’ayer,  cross-ply  (C°/90°)  rect- 
angular plate  uncer  sinusoidal  loading. 


CONCLUDING  REMARKS 

A theory  of  plate  bending  for  laminated,  anistropic  bimodulus  materials 
were  presented.  This  theory  is  based  on  the  fiber-governed  symmetric  compliance 
model  for  bimodulus  materials. 

Closed-form  and  finite-element  solutions  have  been  presented  for  selected 
problems  involving  different  plate  geometries,  lamination  arrangements, 
boundary  conditions,  both  without  and  with  thickness-shear  deformation. 

It  was  shown  that  even  *or  sing'e-layer  bimodulus  p’ates,  the  neutral-sur- 
face location  ^ay  vary  considerably  from  the  geometric  midplane,  defending  upon 
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Fig.  6 Transverse  deflection  vs.  plate  aspect  ratio, 
and  side-to-thickness  ratio  for  single-layer 
and  two-layer  cross-ply  plates  under  sinusoidal 
loading. 
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tho  degree  of  bi  ..odulari  ty.  Also,  the  plate  deflection  is  signi ficantly 
affected  by  tne  biniodulus  action. 
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